####Replication for Nelson and Driscoll, "Accountability for Court Packing"

####REPLICATION FOR FIGURES IN MAIN TEXT AND APPENDIX E

######CCES#######
###Figure 1
##Index DV
setEPS()
postscript("ccesJLC.eps",width=17,height=5,paper='special')
#pdf("ccesJLC.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#0 |   .3693525   .0108469       .348067    .3906379
#1 |   .4956801   .0106047      .4748699    .5164903
plot(NA,ylab="Proposer Support", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.4956801,x=1,pch=15,col="black",cex=1.5)
points(y=.3693525,x=2,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4748699,y1=.5164903,lwd=2, col="black")
segments(x0=2,x1=2,y0=.348067,y1=.3906379,lwd=2, col="black")

#0 |   .3175458   .0115459       .294883    .3402085
#1 |   .5596751   .0105956      .5388775    .5804727
plot(NA,ylab="Proposer Support", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.3175458,x=1,pch=19,col="black",cex=1.5)
points(y=.5596751,x=2,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.294883,y1=.3402085,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5388775,y1=.5804727,lwd=2, col="black")


#1 |    .617887   .0141741      .5900655    .6457085
#2 |   .5026216   .0146771      .4738127    .5314305
#3 |   .3953142   .0154932      .3649034     .425725
#4 |   .2470682   .0155415      .2165625    .2775738
plot(NA,ylab="Proposer Support", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.3953142,x=.9,pch=15,col="black",cex=1.5)
points(y=.2470682,x=1.1,pch=17,col="black",cex=1.5)
points(y=.617887,x=1.9,pch=15,col="black",cex=1.5)
points(y=.5026216,x=2.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.3649034,y1=.425725,lwd=2, col="black")
segments(x0=1.1,x1=1.1,y0=.2165625,y1=.2775738,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.5900655,y1=.6457085,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.4738127,y1=.5314305,lwd=2, col="black")
legend("topleft", legend=c("Bureaucratic", "Politicized"), pch=c(15,17),
       col=c("black", "black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,ylab="Proposer Support", col="black",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy", xlim=c(0,1),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19)
#axis(1, at=c(0,.25,.5,.75,1),cex.axis=1.5,las=1)
lines(x=c(0,.25,.5,.75,1),y=c(0.3555196,
                              0.3371489,
                              0.3187782,
                              0.3004074,
                              0.2820367),lwd=2,lty=2,col="gray60")
lines(x=c(0,.25,.5,.75,1),y=c(0.5491991,
                              0.5540439,
                              0.5588887,
                              0.5637336,
                              0.5685784),lwd=2)
legend("bottomleft", legend=c("Copartisan", "Outpartisan"), lty=c(1,2),
       col=c("black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()

###Figure E1
##Proposal DV
pdf("ccesProposal.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#0 |   .3556202   .0124392      .3312102    .3800301
#1 |   .5217391   .0121331      .4979299    .5455484
plot(NA,ylab="Proposal Support", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.5217391,x=1,pch=15,col="black",cex=1.5)
points(y=.3556202,x=2,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4979299,y1=.5455484,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3312102,y1=.3800301,lwd=2, col="black")

#0 |   .3368298   .0140418      .3092682    .3643915
#1 |   .5523691   .0127989      .5272471    .5774911
plot(NA,ylab="Proposal Support", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.3368298,x=1,pch=19,col="black",cex=1.5)
points(y=.5523691,x=2,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.3092682,y1=.3643915,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5272471,y1=.5774911,lwd=2, col="black")

#1 |      .6275   .0168117      .5945016    .6604984
#2 |   .4776119    .017826      .4426225    .5126014
#3 |   .4411765   .0190599      .4037652    .4785878
#4 |   .2422222   .0183246      .2062541    .2781903

plot(NA,ylab="Proposal Support", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.4411765,x=.9,pch=15,col="black",cex=1.5)
points(y=.2422222,x=1.1,pch=17,col="black",cex=1.5)
points(y=.6275,x=1.9,pch=15,col="black",cex=1.5)
points(y=.4776119,x=2.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.4037652,y1=.4785878,lwd=2, col="black")
segments(x0=1.1,x1=1.1,y0=.2062541,y1=.2781903,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.5945016,y1=.6604984,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.4426225,y1=.5126014,lwd=2, col="black")
legend("topleft", legend=c("Bureaucratic", "Politicized"), pch=c(15,17),
       col=c("black", "black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,ylab="Proposal Support", col="black",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy", xlim=c(0,1),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19)
#axis(1, at=c(0,.25,.5,.75,1),cex.axis=1.5,las=1)
lines(x=c(0,.25,.5,.75,1),y=c(0.3765918,
                              0.3574053,
                              0.3382189,
                              0.3190324,
                              0.299846),lwd=2,lty=2,col="gray60")
lines(x=c(0,.25,.5,.75,1),y=c(0.5478161,
                              0.5478571,
                              0.547898,
                              0.547939,
                              0.5479799),lwd=2)
legend("bottomleft", legend=c("Copartisan", "Outpartisan"), lty=c(1,2),
       col=c("black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()

##Approval DV
pdf("ccesApproval.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#0 |   .3902913   .0109175      .3688672    .4117153
#1 |   .4796451   .0108571      .4583397    .5009505
plot(NA,ylab="Job Performance", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.4796451,x=1,pch=15,col="black",cex=1.5)
points(y=.3902913,x=2,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4583397,y1=.5009505,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3688672,y1=.4117153,lwd=2, col="black")

#0 |   .3407494   .0120804      .3170375    .3644613
#1 |   .5401003    .011041      .5184286    .5617719
plot(NA,ylab="Job Performance", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.3407494,x=1,pch=19,col="black",cex=1.5)
points(y=.5401003,x=2,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.3170375,y1=.3644613,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5184286,y1=.5617719,lwd=2, col="black")


#1 |   .5824873   .0152435      .5525667     .612408
#2 |   .4987624   .0154463      .4684437    .5290811
#3 |   .3977833   .0167039      .3649961    .4305704
#4 |   .2890625   .0166487      .2563836    .3217414

plot(NA,ylab="Job Performance", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.3977833,x=.9,pch=15,col="black",cex=1.5)
points(y=.2890625,x=1.1,pch=17,col="black",cex=1.5)
points(y=.5824873,x=1.9,pch=15,col="black",cex=1.5)
points(y=.4987624,x=2.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.3649961,y1=.4305704,lwd=2, col="black")
segments(x0=1.1,x1=1.1,y0=.2563836,y1=.3217414,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.5525667,y1=.612408,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.4684437,y1=.5290811,lwd=2, col="black")
legend("topleft", legend=c("Bureaucratic", "Politicized"), pch=c(15,17),
       col=c("black", "black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,ylab="Job Performance", col="black",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy", xlim=c(0,1),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19)
#axis(1, at=c(0,.25,.5,.75,1),cex.axis=1.5,las=1)
lines(x=c(0,.25,.5,.75,1),y=c(0.3772684,
                              0.3601127,
                              0.3429569,
                              0.3258011,
                              0.3086453),lwd=2,lty=2,col="gray60")
lines(x=c(0,.25,.5,.75,1),y=c(0.5380046,
                              0.5397316,
                              0.5414586,
                              0.5431856,
                              0.5449126),lwd=2)
legend("bottomleft", legend=c("Copartisan", "Outpartisan"), lty=c(1,2),
       col=c("black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()

##Vote DV
pdf("ccesVote.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#0 |   .3597679   .0143041      .3316984    .3878373
#1 |   .4782609   .0146583      .4494964    .5070254
plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.4782609,x=1,pch=15,col="black",cex=1.5)
points(y=.3597679,x=2,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4494964,y1=.5070254,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3316984,y1=.3878373,lwd=2, col="black")

#0 |   .2517483   .0139656      .2243363    .2791602
#1 |   .6026119   .0146462      .5738639    .6313599
plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.2517483,x=1,pch=19,col="black",cex=1.5)
points(y=.6026119,x=2,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.2243363,y1=.2791602,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5738639,y1=.6313599,lwd=2, col="black")

#1 |       .655   .0204042      .6149501    .6950499
#2 |   .5507426   .0204123      .5106767    .5908085
#3 |   .3161765   .0198739      .2771674    .3551855
#4 |   .1933333    .018814      .1564046     .230262
plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,2.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2),labels=c("Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.3161765,x=.9,pch=15,col="black",cex=1.5)
points(y=.1933333,x=1.1,pch=17,col="black",cex=1.5)
points(y=.655,x=1.9,pch=15,col="black",cex=1.5)
points(y=.5507426,x=2.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.2771674,y1=.3551855,lwd=2, col="black")
segments(x0=1.1,x1=1.1,y0=.1564046,y1=.230262,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.6149501,y1=.6950499,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.5106767,y1=.5908085,lwd=2, col="black")
legend("topleft", legend=c("Bureaucratic", "Politicized"), pch=c(15,17),
       col=c("black", "black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy", xlim=c(0,1),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19)
#axis(1, at=c(0,.25,.5,.75,1),cex.axis=1.5,las=1)
lines(x=c(0,.25,.5,.75,1),y=c(0.292591,
                              0.2719987,
                              0.2514064,
                              0.2308141,
                              0.2102218),lwd=2,lty=2,col="gray60")
lines(x=c(0,.25,.5,.75,1),y=c(0.5650563,
                              0.5835552,
                              0.602054,
                              0.6205529,
                              0.6390518),lwd=2)
legend("bottomleft", legend=c("Copartisan", "Outpartisan"), lty=c(1,2),
       col=c("black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()


###MTURK1###
###Figure 2
##Index DV
setEPS()
postscript("Turk1JLC.eps",width=17,height=5)
#pdf("Turk1JLC.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#Control |   .4680874   .0083736      .4516673    .4845074
#Bureaucratic |   .5893621   .0079066      .5738578    .6048664
#Politicized |   .3770168    .009137      .3590998    .3949337
plot(NA,ylab="Proposer Support", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("No Rationale","Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.4680874,x=1,pch=19,col="black",cex=1.5)
points(y=.5893621,x=2,pch=15,col="black",cex=1.5)
points(y=.3770168,x=3,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4516673,y1=.4845074,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5738578,y1=.6048664,lwd=2, col="black")
segments(x0=3,x1=3,y0=.3590998,y1=.3949337,lwd=2, col="black")

#Control |   .4754161   .0110235      .4537965    .4970358
#Outpartisan |   .4474259   .0104972      .4268385    .4680134
#Copartisan |   .5406215   .0096513      .5216931      .55955
plot(NA,ylab="Proposer Support", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.4754161,x=1,pch=19,col="black",cex=1.5)
points(y=.4474259,x=2,pch=19,col="black",cex=1.5)
points(y=.5406215,x=3,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4537965,y1=.4970358,lwd=2, col="black")
segments(x0=2,x1=2,y0=.4268385,y1=.4680134,lwd=2, col="black")
segments(x0=3,x1=3,y0=.5216931,y1=.55955,lwd=2, col="black")


plot(NA,ylab="Proposer Support", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
#CC |   .4653232   .0176168      .4307726    .4998738
#CCo |   .5400947   .0152816       .510124    .5700654
#COth |   .4202881   .0169162      .3871115    .4534647
#BC |   .5833268   .0169206      .5501417    .6165119
#BCo |   .6420218   .0147465      .6131006     .670943
#BOth |   .5608867   .0154442      .5305971    .5911762
#PC |   .3799707   .0196882      .3413577    .4185838
#PCO |   .4361732   .0174118      .4020247    .4703218
#POth |   .3497299   .0194568      .3115707     .387889
points(y=.4653232,x=.9,pch=19,col="black",cex=1.5)
points(y=.4202881,x=1.9,pch=19,col="black",cex=1.5)
points(y=.5400947,x=2.9,pch=19,col="black",cex=1.5)

points(y=.5833268,x=1,pch=15,col="black",cex=1.5)
points(y=.5608867,x=2,pch=15,col="black",cex=1.5)
points(y=.6420218,x=3,pch=15,col="black",cex=1.5)

points(y=.3799707,x=1.1,pch=17,col="black",cex=1.5)
points(y=.3497299,x=2.1,pch=17,col="black",cex=1.5)
points(y=.4361732,x=3.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.4307726,y1=.4998738,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.3871115,y1=.4534647,lwd=2, col="black")
segments(x0=2.9,x1=2.9,y0=.510124,y1=.5700654,lwd=2, col="black")

segments(x0=1,x1=1,y0=.5501417,y1=.6165119,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5305971,y1=.5911762,lwd=2, col="black")
segments(x0=3,x1=3,y0=.6131006,y1=.670943,lwd=2, col="black")

segments(x0=1.1,x1=1.1,y0=.3413577,y1=.4185838,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.3115707,y1=.387889,lwd=2, col="black")
segments(x0=3.1,x1=3.1,y0=.4020247,y1=.4703218,lwd=2, col="black")
legend("bottomright", legend=c("No Rationale","Bureaucratic", "Politicized"), pch=c(19,15,17),
       col=c("black", "black","black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))


plot(NA,xlim=c(0,1),ylim=c(0,.8),ylab="Proposer Support", col="#0072ce",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy",
     cex.axis=1.5,cex.lab=2)
#axis(1, at=c(0,.2,.4,.6,.8,1),labels=c("0.00","0.20","0.40","0.60","0.80","1.00"),cex.axis=1.5)
lines(x=seq(0,1,by=.05),y=c(0.5462058,
                            0.5394465,
                            0.5326872,
                            0.5259279,
                            0.5191686,
                            0.5124093,
                            0.50565,
                            0.4988907,
                            0.4921314,
                            0.4853721,
                            0.4786128,
                            0.4718534,
                            0.4650941,
                            0.4583348,
                            0.4515755,
                            0.4448162,
                            0.4380569,
                            0.4312976,
                            0.4245383,
                            0.417779,
                            0.4110197),lwd=2,lty=3,col="gray30")
lines(x=seq(0,1,by=.05),y=c(0.5445685,
                            0.543862,
                            0.5431555,
                            0.542449,
                            0.5417425,
                            0.541036,
                            0.5403295,
                            0.5396229,
                            0.5389164,
                            0.5382099,
                            0.5375034,
                            0.5367969,
                            0.5360904,
                            0.5353839,
                            0.5346774,
                            0.5339709,
                            0.5332644,
                            0.5325579,
                            0.5318514,
                            0.5311449,
                            0.5304383),lwd=2)
lines(x=seq(0,1,by=.05),y=c(0.5381136,
                            0.5293171,
                            0.5205206,
                            0.5117242,
                            0.5029277,
                            0.4941312,
                            0.4853348,
                            0.4765383,
                            0.4677419,
                            0.4589454,
                            0.4501489,
                            0.4413525,
                            0.432556,
                            0.4237595,
                            0.4149631,
                            0.4061666,
                            0.3973701,
                            0.3885737,
                            0.3797772,
                            0.3709807,
                            0.3621843),lwd=2,lty=2,col="gray60")
legend("bottomleft", legend=c("Nonpartisan","Copartisan", "Outpartisan"), lty=c(3,1,2),
       col=c("gray30","black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()

###Figure E2
##Proposal
pdf("Turk1Proposal.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#Control |   .5021135   .0094526      .4835777    .5206494
#Bureaucratic |   .6438107   .0085075       .627128    .6604934
#Politicized |   .3864198   .0102745      .3662722    .4065673
plot(NA,ylab="Proposal Support", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("No Rationale","Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.5021135,x=1,pch=19,col="black",cex=1.5)
points(y=.6438107,x=2,pch=15,col="black",cex=1.5)
points(y=.3864198,x=3,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4835777,y1=.5206494,lwd=2, col="black")
segments(x0=2,x1=2,y0=.627128,y1=.6604934,lwd=2, col="black")
segments(x0=3,x1=3,y0=.3662722,y1=.4065673,lwd=2, col="black")

#Control |   .5145985   .0121853      .4907004    .5384966
#Copartisan |   .5659026   .0106898      .5449376    .5868676
#Outpartisan |    .478482   .0119165      .4551112    .5018528

plot(NA,ylab="Proposal Support", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.5145985,x=1,pch=19,col="black",cex=1.5)
points(y=.478482,x=2,pch=19,col="black",cex=1.5)
points(y=.5659026,x=3,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4907004,y1=.5384966,lwd=2, col="black")
segments(x0=2,x1=2,y0=.4551112,y1=.5018528,lwd=2, col="black")
segments(x0=3,x1=3,y0=.5449376,y1=.5868676,lwd=2, col="black")


plot(NA,ylab="Proposal Support", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
#CC |   .5105263   .0198181      .4716585    .5493941
#CCo |   .5588235   .0170001      .5254826    .5921645
#COth |   .4622222   .0193797      .4242143    .5002302
#BC |    .630618    .017727      .5958514    .6653846
#BCo |   .6987179   .0152979      .6687153    .7287206
#BOth |   .6131222   .0169974      .5797864     .646458
#PC |   .4041667    .022221      .3605863     .447747
#PCO |   .4358407   .0191098       .398362    .4733194
#POth |   .3432642   .0216202      .3008621    .3856664

points(y=.5105263,x=.9,pch=19,col="black",cex=1.5)
points(y=.4622222,x=1.9,pch=19,col="black",cex=1.5)
points(y=.5588235,x=2.9,pch=19,col="black",cex=1.5)

points(y=.630618,x=1,pch=15,col="black",cex=1.5)
points(y=.6131222,x=2,pch=15,col="black",cex=1.5)
points(y=.6987179,x=3,pch=15,col="black",cex=1.5)

points(y=.4041667,x=1.1,pch=17,col="black",cex=1.5)
points(y=.3432642,x=2.1,pch=17,col="black",cex=1.5)
points(y=.4358407,x=3.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.4716585,y1=.5493941,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.4242143,y1=.5002302,lwd=2, col="black")
segments(x0=2.9,x1=2.9,y0=.5254826,y1=.5921645,lwd=2, col="black")

segments(x0=1,x1=1,y0=.5958514,y1=.6653846,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5797864,y1=.646458,lwd=2, col="black")
segments(x0=3,x1=3,y0=.6687153,y1=.7287206,lwd=2, col="black")

segments(x0=1.1,x1=1.1,y0=.3605863,y1=.447747,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.3008621,y1=.3856664,lwd=2, col="black")
segments(x0=3.1,x1=3.1,y0=.398362,y1=.4733194,lwd=2, col="black")
legend("bottomright", legend=c("No Rationale","Bureaucratic", "Politicized"), pch=c(19,15,17),
       col=c("black", "black","black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,xlim=c(0,1),ylim=c(0,.8),ylab="Proposal Support", col="#0072ce",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy",
     cex.axis=1.5,cex.lab=2)
#axis(1, at=c(0,.2,.4,.6,.8,1),labels=c("0.00","0.20","0.40","0.60","0.80","1.00"),cex.axis=1.5)
lines(x=seq(0,1,by=.05),y=c(0.5884137,
                            0.5811564,
                            0.5738991,
                            0.5666417,
                            0.5593844,
                            0.5521271,
                            0.5448698,
                            0.5376125,
                            0.5303551,
                            0.5230978,
                            0.5158405,
                            0.5085832,
                            0.5013258,
                            0.4940685,
                            0.4868112,
                            0.4795539,
                            0.4722966,
                            0.4650392,
                            0.4577819,
                            0.4505246,
                            0.4432673),lwd=2,lty=3,col="gray30")
lines(x=seq(0,1,by=.05),y=c(0.5546377,
                            0.5552351,
                            0.5558326,
                            0.55643,
                            0.5570275,
                            0.5576249,
                            0.5582224,
                            0.5588198,
                            0.5594172,
                            0.5600147,
                            0.5606121,
                            0.5612096,
                            0.561807,
                            0.5624045,
                            0.5630019,
                            0.5635993,
                            0.5641968,
                            0.5647942,
                            0.5653917,
                            0.5659891,
                            0.5665866),lwd=2)
lines(x=seq(0,1,by=.05),y=c(0.5519592,
                            0.5443793,
                            0.5367994,
                            0.5292196,
                            0.5216397,
                            0.5140599,
                            0.50648,
                            0.4989002,
                            0.4913203,
                            0.4837404,
                            0.4761606,
                            0.4685807,
                            0.4610009,
                            0.453421,
                            0.4458411,
                            0.4382613,
                            0.4306814,
                            0.4231016,
                            0.4155217,
                            0.4079419,
                            0.400362),lwd=2,lty=2,col="gray60")
legend("bottomleft", legend=c("Nonpartisan","Copartisan", "Outpartisan"), lty=c(3,1,2),
       col=c("gray30","black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()


###Incumbent
pdf("Turk1Approval.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#Control |    .509047   .0083346      .4927035    .5253906
#Bureaucratic |   .5932927   .0076923      .5782087    .6083767
#Politicized |   .4371147   .0092673      .4189421    .4552873
plot(NA,ylab="Job Performance", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("No Rationale","Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.509047,x=1,pch=19,col="black",cex=1.5)
points(y=.5932927,x=2,pch=15,col="black",cex=1.5)
points(y=.4371147,x=3,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4927035,y1=.5253906,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5782087,y1=.6083767,lwd=2, col="black")
segments(x0=3,x1=3,y0=.4189421,y1=.4552873,lwd=2, col="black")

#Control |   .5142202   .0109351      .4927739    .5356665
#Copartisan |      .5675   .0091397       .549575     .585425
#Outpartisan |   .4816119   .0104149      .4611859    .5020379

plot(NA,ylab="Job Performance", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.5142202,x=1,pch=19,col="black",cex=1.5)
points(y=.4816119,x=2,pch=19,col="black",cex=1.5)
points(y=.5675,x=3,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4927739,y1=.5356665,lwd=2, col="black")
segments(x0=2,x1=2,y0=.4611859,y1=.5020379,lwd=2, col="black")
segments(x0=3,x1=3,y0=.549575,y1=.585425,lwd=2, col="black")


plot(NA,ylab="Job Performance", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
#CC |   .5026455   .0178912      .4675567    .5377343
#CCo |     .58159   .0149181      .5523323    .6108476
#COth |   .4535398   .0175361      .4191475    .4879321
#BC |   .6036932   .0159087      .5724926    .6348938
#BCo |   .6260684   .0148547      .5969349    .6552018
#BOth |   .5693182   .0152844      .5393421    .5992942
#PC |   .4388889   .0207164      .3982595    .4795183
#PCO |   .4922907   .0165484      .4598356    .5247459
#POth |   .4145078   .0198422      .3755928    .4534227

points(y=.5026455,x=.9,pch=19,col="black",cex=1.5)
points(y=.4535398,x=1.9,pch=19,col="black",cex=1.5)
points(y=.58159,x=2.9,pch=19,col="black",cex=1.5)

points(y=.6036932,x=1,pch=15,col="black",cex=1.5)
points(y=.5693182,x=2,pch=15,col="black",cex=1.5)
points(y=.6260684,x=3,pch=15,col="black",cex=1.5)

points(y=.4388889,x=1.1,pch=17,col="black",cex=1.5)
points(y=.4145078,x=2.1,pch=17,col="black",cex=1.5)
points(y=.4922907,x=3.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.4675567,y1=.5377343,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.4191475,y1=.4879321,lwd=2, col="black")
segments(x0=2.9,x1=2.9,y0=.5523323,y1=.6108476,lwd=2, col="black")

segments(x0=1,x1=1,y0=.5724926,y1=.6348938,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5393421,y1=.5992942,lwd=2, col="black")
segments(x0=3,x1=3,y0=.5969349,y1=.6552018,lwd=2, col="black")

segments(x0=1.1,x1=1.1,y0=.3982595,y1=.4795183,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.3755928,y1=.4534227,lwd=2, col="black")
segments(x0=3.1,x1=3.1,y0=.4598356,y1=.5247459,lwd=2, col="black")
legend("bottomright", legend=c("No Rationale","Bureaucratic", "Politicized"), pch=c(19,15,17),
       col=c("black", "black","black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,xlim=c(0,1),ylim=c(0,.8),ylab="Job Performance", col="#0072ce",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy",
     cex.axis=1.5,cex.lab=2)
#axis(1, at=c(0,.2,.4,.6,.8,1),labels=c("0.00","0.20","0.40","0.60","0.80","1.00"),cex.axis=1.5)
lines(x=seq(0,1,by=.05),y=c(0.5834994,
                            0.5767621,
                            0.5700248,
                            0.5632875,
                            0.5565503,
                            0.549813,
                            0.5430757,
                            0.5363384,
                            0.5296012,
                            0.5228639,
                            0.5161266,
                            0.5093893,
                            0.5026521,
                            0.4959148,
                            0.4891775,
                            0.4824403,
                            0.475703,
                            0.4689657,
                            0.4622284,
                            0.4554912,
                            0.4487539),lwd=2,lty=3,col="gray30")
lines(x=seq(0,1,by=.05),y=c(0.5865023,
                            0.5844307,
                            0.582359,
                            0.5802874,
                            0.5782158,
                            0.5761441,
                            0.5740725,
                            0.5720009,
                            0.5699293,
                            0.5678576,
                            0.565786,
                            0.5637144,
                            0.5616427,
                            0.5595711,
                            0.5574995,
                            0.5554279,
                            0.5533562,
                            0.5512846,
                            0.549213,
                            0.5471413,
                            0.5450697),lwd=2)
lines(x=seq(0,1,by=.05),y=c(0.5684044,
                            0.5600279,
                            0.5516514,
                            0.543275,
                            0.5348985,
                            0.526522,
                            0.5181456,
                            0.5097691,
                            0.5013926,
                            0.4930162,
                            0.4846397,
                            0.4762632,
                            0.4678868,
                            0.4595103,
                            0.4511339,
                            0.4427574,
                            0.4343809,
                            0.4260045,
                            0.417628,
                            0.4092515,
                            0.4008751),lwd=2,lty=2,col="gray60")
legend("bottomleft", legend=c("Nonpartisan","Copartisan", "Outpartisan"), lty=c(3,1,2),
       col=c("gray30","black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))


dev.off()


###Vote
pdf("Turk1Vote.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#Control |   .2373494   .0147768      .2083732    .2663256
#Bureaucratic |   .4058323   .0171274      .3722467     .439418
#Politicized |   .1742892   .0133458      .1481191    .2004593
plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("No Rationale","Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.2373494,x=1,pch=19,col="black",cex=1.5)
points(y=.4058323,x=2,pch=15,col="black",cex=1.5)
points(y=.1742892,x=3,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.2083732,y1=.2663256,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3722467,y1=.439418,lwd=2, col="black")
segments(x0=3,x1=3,y0=.1481191,y1=.2004593,lwd=2, col="black")

#Control |   .2413163   .0183116       .205403    .2772295
#Copartisan |   .3796562   .0183821      .3436047    .4157076
#Outpartisan |    .253125   .0172005       .219391     .286859

plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.2413163,x=1,pch=19,col="black",cex=1.5)
points(y=.253125,x=2,pch=19,col="black",cex=1.5)
points(y=.3796562,x=3,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.205403,y1=.2772295,lwd=2, col="black")
segments(x0=2,x1=2,y0=.219391,y1=.286859,lwd=2, col="black")
segments(x0=3,x1=3,y0=.3436047,y1=.4157076,lwd=2, col="black")

plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
#CC |         .2   .0290957      .1429368    .2570632
#CCo |   .3556485   .0310301      .2947915    .4165055
#COth |   .1946903   .0263975      .1429189    .2464616
#BC |   .3820225   .0365211      .3103964    .4536486
#BCo |   .5107296    .032819      .4463641    .5750951
#BOth |   .3755656   .0326494      .3115329    .4395984
#PC |   .1452514   .0264101      .0934554    .1970474
#PCO |   .2699115   .0295942      .2118706    .3279524
#POth |   .1813472    .027807      .1268113     .235883

points(y=.2,x=.9,pch=19,col="black",cex=1.5)
points(y=.1946903,x=1.9,pch=19,col="black",cex=1.5)
points(y=.3556485,x=2.9,pch=19,col="black",cex=1.5)

points(y=.3820225,x=1,pch=15,col="black",cex=1.5)
points(y=.3755656,x=2,pch=15,col="black",cex=1.5)
points(y=.5107296,x=3,pch=15,col="black",cex=1.5)

points(y=.1452514,x=1.1,pch=17,col="black",cex=1.5)
points(y=.1813472,x=2.1,pch=17,col="black",cex=1.5)
points(y=.2699115,x=3.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.1429368,y1=.2570632,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.1429189,y1=.2464616,lwd=2, col="black")
segments(x0=2.9,x1=2.9,y0=.2947915,y1=.4165055,lwd=2, col="black")

segments(x0=1,x1=1,y0=.3103964,y1=.4536486,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3115329,y1=.4395984,lwd=2, col="black")
segments(x0=3,x1=3,y0=.4463641,y1=.5750951,lwd=2, col="black")

segments(x0=1.1,x1=1.1,y0=.0934554,y1=.1970474,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.1268113,y1=.235883,lwd=2, col="black")
segments(x0=3.1,x1=3.1,y0=.2118706,y1=.3279524,lwd=2, col="black")
legend("bottomright", legend=c("No Rationale","Bureaucratic", "Politicized"), pch=c(19,15,17),
       col=c("black", "black","black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,xlim=c(0,1),ylim=c(0,.8),ylab="Vote for Incumbent?", col="#0072ce",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy",
     cex.axis=1.5,cex.lab=2)
#axis(1, at=c(0,.2,.4,.6,.8,1),labels=c("0.00","0.20","0.40","0.60","0.80","1.00"),cex.axis=1.5)
lines(x=seq(0,1,by=.05),y=c(0.3155899,
                            0.3096209,
                            0.303652,
                            0.2976831,
                            0.2917141,
                            0.2857452,
                            0.2797762,
                            0.2738073,
                            0.2678383,
                            0.2618694,
                            0.2559005,
                            0.2499315,
                            0.2439626,
                            0.2379936,
                            0.2320247,
                            0.2260557,
                            0.2200868,
                            0.2141178,
                            0.2081489,
                            0.20218,
                            0.196211),lwd=2,lty=3,col="gray30")
lines(x=seq(0,1,by=.05),y=c(0.3763368,
                            0.3766884,
                            0.3770399,
                            0.3773915,
                            0.3777431,
                            0.3780947,
                            0.3784462,
                            0.3787978,
                            0.3791494,
                            0.3795009,
                            0.3798525,
                            0.3802041,
                            0.3805556,
                            0.3809072,
                            0.3812588,
                            0.3816103,
                            0.3819619,
                            0.3823135,
                            0.3826651,
                            0.3830166,
                            0.3833682),lwd=2)
lines(x=seq(0,1,by=.05),y=c(0.4159744,
                            0.4014993,
                            0.3870242,
                            0.3725491,
                            0.358074,
                            0.3435989,
                            0.3291238,
                            0.3146487,
                            0.3001736,
                            0.2856985,
                            0.2712234,
                            0.2567483,
                            0.2422732,
                            0.2277981,
                            0.213323,
                            0.1988478,
                            0.1843727,
                            0.1698976,
                            0.1554225,
                            0.1409474,
                            0.1264723),lwd=2,lty=2,col="gray60")
legend("bottomleft", legend=c("Nonpartisan","Copartisan", "Outpartisan"), lty=c(3,1,2),
       col=c("gray30","black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()


###MTURK2###
###Figure 3
##Index
setEPS()
postscript("Turk2JLC.eps",width=17,height=5)
#pdf("Turk2JLC.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#Control |   .4680259    .015897       .436828    .4992238
#Bureaucratic |   .5281718   .0155679      .4976197    .5587239
#Politicized |   .4232349    .016219       .391405    .4550648
plot(NA,ylab="Proposer Support", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("No Rationale","Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.4680259,x=1,pch=19,col="black",cex=1.5)
points(y=.5281718,x=2,pch=15,col="black",cex=1.5)
points(y=.4232349,x=3,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.436828,y1=.4992238,lwd=2, col="black")
segments(x0=2,x1=2,y0=.4976197,y1=.5587239,lwd=2, col="black")
segments(x0=3,x1=3,y0=.391405,y1=.4550648,lwd=2, col="black")


#Control |   .4960948   .0166877      .4633412    .5288484
#Outpartisan |   .3670287   .0169026      .3338533    .4002041
#Copartisan |   .6004582   .0139578      .5730625    .6278538
plot(NA,ylab="Proposer Support", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.4960948,x=1,pch=19,col="black",cex=1.5)
points(y=.3670287,x=2,pch=19,col="black",cex=1.5)
points(y=.6004582,x=3,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4633412,y1=.5288484,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3338533,y1=.4002041,lwd=2, col="black")
segments(x0=3,x1=3,y0=.5730625,y1=.6278538,lwd=2, col="black")

plot(NA,ylab="Proposer Support", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
#CC |   .4706079   .0268968      .4178162    .5233995
#CCo |   .5848871   .0254639       .534908    .6348663
#COth |     .37842   .0312484      .3170872    .4397528
#BC |   .5795434   .0268889      .5267674    .6323194
#BCo |   .6281156    .021609      .5857026    .6705286
#BOth |   .4238209   .0304002       .364153    .4834889
#PC |   .4429506   .0311112      .3818872    .5040139
#PCO |   .5889305   .0248481      .5401599    .6377011
#POth |   .3148311    .026148      .2635093    .3661529
points(y=.4706079,x=.9,pch=19,col="black",cex=1.5)
points(y=.37842,x=1.9,pch=19,col="black",cex=1.5)
points(y=.5848871,x=2.9,pch=19,col="black",cex=1.5)
points(y=.5795434,x=1,pch=15,col="black",cex=1.5)
points(y=.4238209,x=2,pch=15,col="black",cex=1.5)
points(y=.6281156,x=3,pch=15,col="black",cex=1.5)
points(y=.4429506,x=1.1,pch=17,col="black",cex=1.5)
points(y=.3148311,x=2.1,pch=17,col="black",cex=1.5)
points(y=.5889305,x=3.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.4178162,y1=.5233995,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.3170872,y1=.4397528,lwd=2, col="black")
segments(x0=2.9,x1=2.9,y0=.534908,y1=.6348663,lwd=2, col="black")
segments(x0=1,x1=1,y0=.5267674,y1=.6323194,lwd=2, col="black")
segments(x0=2,x1=2,y0=.364153,y1=.4834889,lwd=2, col="black")
segments(x0=3,x1=3,y0=.5857026,y1=.6705286,lwd=2, col="black")
segments(x0=1.1,x1=1.1,y0=.3818872,y1=.5040139,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.2635093,y1=.3661529,lwd=2, col="black")
segments(x0=3.1,x1=3.1,y0=.5401599,y1=.6377011,lwd=2, col="black")
legend("bottomright", legend=c("No Rationale","Bureaucratic", "Politicized"), pch=c(19,15,17),
       col=c("black", "black","black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,xlim=c(0,1),ylim=c(0,.8),ylab="Proposer Support", col="#0072ce",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy",
     cex.axis=1.5,cex.lab=2)
#axis(1, at=c(0,.2,.4,.6,.8,1),labels=c("0.00","0.20","0.40","0.60","0.80","1.00"),cex.axis=1.5)
lines(x=seq(0,1,by=.05),y=c(0.5698168,
                            0.5618556,
                            0.5538944,
                            0.5459332,
                            0.537972,
                            0.5300108,
                            0.5220496,
                            0.5140884,
                            0.5061272,
                            0.498166,
                            0.4902048,
                            0.4822436,
                            0.4742823,
                            0.4663211,
                            0.4583599,
                            0.4503987,
                            0.4424375,
                            0.4344763,
                            0.4265151,
                            0.4185539,
                            0.4105927),lwd=2,lty=3,col="gray30")
lines(x=seq(0,1,by=.05),y=c(0.5730654,
                            0.5757573,
                            0.5784491,
                            0.5811409,
                            0.5838328,
                            0.5865246,
                            0.5892164,
                            0.5919083,
                            0.5946001,
                            0.597292,
                            0.5999838,
                            0.6026756,
                            0.6053675,
                            0.6080593,
                            0.6107511,
                            0.613443,
                            0.6161348,
                            0.6188267,
                            0.6215185,
                            0.6242103,
                            0.6269022),lwd=2)
lines(x=seq(0,1,by=.05),y=c(0.547114,
                            0.5318749,
                            0.5166358,
                            0.5013966,
                            0.4861575,
                            0.4709184,
                            0.4556793,
                            0.4404402,
                            0.4252011,
                            0.4099619,
                            0.3947228,
                            0.3794837,
                            0.3642446,
                            0.3490055,
                            0.3337664,
                            0.3185272,
                            0.3032881,
                            0.288049,
                            0.2728099,
                            0.2575708,
                            0.2423317),lwd=2,lty=2,col="gray60")
legend("bottomleft", legend=c("Nonpartisan","Copartisan", "Outpartisan"), lty=c(3,1,2),
       col=c("gray30","black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()

###Proposal
pdf("Turk2Proposal.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#Control |   .4782245   .0144712      .4498334    .5066156
#Bureaucratic |    .539644   .0133884      .5133774    .5659107
#Politicized |   .3995215   .0145151      .3710443    .4279988
plot(NA,ylab="Proposal Support", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("No Rationale","Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.4782245,x=1,pch=19,col="black",cex=1.5)
points(y=.539644,x=2,pch=15,col="black",cex=1.5)
points(y=.3995215,x=3,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4498334,y1=.5066156,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5133774,y1=.5659107,lwd=2, col="black")
segments(x0=3,x1=3,y0=.3710443,y1=.4279988,lwd=2, col="black")

#Control |   .5008818   .0140311      .4733515    .5284122
#Copartisan |   .5833333   .0137424      .5563693    .6102974
#Outpartisan |   .3717278    .015407      .3414977    .4019578
plot(NA,ylab="Proposal Support", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.5008818,x=1,pch=19,col="black",cex=1.5)
points(y=.3717278,x=2,pch=19,col="black",cex=1.5)
points(y=.5833333,x=3,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4733515,y1=.5284122,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3414977,y1=.4019578,lwd=2, col="black")
segments(x0=3,x1=3,y0=.5563693,y1=.6102974,lwd=2, col="black")

plot(NA,ylab="Proposal Support", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
#CC |   .4769231   .0241296      .4295784    .5242678
#CCo |   .5820106   .0243248      .5342828    .6297384
#COth |   .3921569   .0284078       .336418    .4478958
#BC |   .5873016   .0210579      .5459838    .6286194
#BCo |   .6318841   .0207389      .5911922    .6725759
#BOth |   .4411028   .0257907      .3904989    .4917067
#PC |   .4371585   .0257603      .3866141    .4877029
#PCO |   .5345345   .0253394      .4848161     .584253
#POth |   .2947846    .024615      .2464875    .3430817

points(y=.4769231,x=.9,pch=19,col="black",cex=1.5)
points(y=.3921569,x=1.9,pch=19,col="black",cex=1.5)
points(y=.5820106,x=2.9,pch=19,col="black",cex=1.5)

points(y=.5873016,x=1,pch=15,col="black",cex=1.5)
points(y=.4411028,x=2,pch=15,col="black",cex=1.5)
points(y=.6318841,x=3,pch=15,col="black",cex=1.5)

points(y=.4371585,x=1.1,pch=17,col="black",cex=1.5)
points(y=.2947846,x=2.1,pch=17,col="black",cex=1.5)
points(y=.5345345,x=3.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.4295784,y1=.5242678,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.336418,y1=.4478958,lwd=2, col="black")
segments(x0=2.9,x1=2.9,y0=.5342828,y1=.6297384,lwd=2, col="black")

segments(x0=1,x1=1,y0=.5459838,y1=.6286194,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3904989,y1=.4917067,lwd=2, col="black")
segments(x0=3,x1=3,y0=.5911922,y1=.6725759,lwd=2, col="black")

segments(x0=1.1,x1=1.1,y0=.3866141,y1=.4877029,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.2464875,y1=.3430817,lwd=2, col="black")
segments(x0=3.1,x1=3.1,y0=.4848161,y1=.584253,lwd=2, col="black")
legend("bottomright", legend=c("No Rationale","Bureaucratic", "Politicized"), pch=c(19,15,17),
       col=c("black", "black","black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,xlim=c(0,1),ylim=c(0,.8),ylab="Proposal Support", col="#0072ce",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy",
     cex.axis=1.5,cex.lab=2)
#axis(1, at=c(0,.2,.4,.6,.8,1),labels=c("0.00","0.20","0.40","0.60","0.80","1.00"),cex.axis=1.5)
lines(x=seq(0,1,by=.05),y=c(0.5624089,
                            0.5567968,
                            0.5511847,
                            0.5455726,
                            0.5399605,
                            0.5343484,
                            0.5287363,
                            0.5231242,
                            0.5175121,
                            0.5119,
                            0.5062879,
                            0.5006758,
                            0.4950637,
                            0.4894516,
                            0.4838395,
                            0.4782274,
                            0.4726153,
                            0.4670032,
                            0.4613911,
                            0.455779,
                            0.4501669),lwd=2,lty=3,col="gray30")
lines(x=seq(0,1,by=.05),y=c(0.5735018,
                            0.5744287,
                            0.5753556,
                            0.5762825,
                            0.5772094,
                            0.5781364,
                            0.5790633,
                            0.5799902,
                            0.5809171,
                            0.5818441,
                            0.582771,
                            0.5836979,
                            0.5846248,
                            0.5855517,
                            0.5864787,
                            0.5874056,
                            0.5883325,
                            0.5892594,
                            0.5901863,
                            0.5911133,
                            0.5920402),lwd=2)
lines(x=seq(0,1,by=.05),y=c(0.5030809,
                            0.4920106,
                            0.4809403,
                            0.46987,
                            0.4587997,
                            0.4477294,
                            0.4366591,
                            0.4255888,
                            0.4145184,
                            0.4034481,
                            0.3923778,
                            0.3813075,
                            0.3702372,
                            0.3591669,
                            0.3480966,
                            0.3370263,
                            0.325956,
                            0.3148856,
                            0.3038153,
                            0.292745,
                            0.2816747),lwd=2,lty=2,col="gray60")
legend("bottomleft", legend=c("Nonpartisan","Copartisan", "Outpartisan"), lty=c(3,1,2),
       col=c("gray30","black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()

###Figure E3
##Incumbent
pdf("Turk2Approval.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#Control |   .4897361   .0146819      .4609247    .5185474
#Bureaucratic |   .5363248     .01522      .5064576     .566192
#Politicized |   .4567164   .0156963      .4259144    .4875184
plot(NA,ylab="Job Performance", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("No Rationale","Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.4897361,x=1,pch=19,col="black",cex=1.5)
points(y=.5363248,x=2,pch=15,col="black",cex=1.5)
points(y=.4567164,x=3,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4609247,y1=.5185474,lwd=2, col="black")
segments(x0=2,x1=2,y0=.5064576,y1=.566192,lwd=2, col="black")
segments(x0=3,x1=3,y0=.4259144,y1=.4875184,lwd=2, col="black")

#Control |    .516317   .0156997      .4855048    .5471293
#Copartisan |   .6079096    .012442       .583491    .6323283
#Outpartisan |   .3981191   .0168739      .3650024    .4312359
plot(NA,ylab="Job Performance", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.516317,x=1,pch=19,col="black",cex=1.5)
points(y=.3981191,x=2,pch=19,col="black",cex=1.5)
points(y=.6079096,x=3,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.4855048,y1=.5471293,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3650024,y1=.4312359,lwd=2, col="black")
segments(x0=3,x1=3,y0=.583491,y1=.6323283,lwd=2, col="black")

plot(NA,ylab="Job Performance", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
#CC |   .4901961   .0248839      .4413588    .5390334
#CCo |   .5958702   .0219547      .5527818    .6389586
#COth |   .4050179   .0296507      .3468254    .4632105
#BC |   .5942029   .0246147      .5458939    .6425119
#BCo |   .6170213    .020762      .5762737    .6577689
#BOth |   .4387755   .0315562      .3768432    .5007078
#PC |   .4673913   .0305059      .4075203    .5272624
#PCO |   .6136364   .0214736      .5714921    .6557806
#POth |   .3619792   .0266688      .3096389    .4143194
points(y=.4901961,x=.9,pch=19,col="black",cex=1.5)
points(y=.4050179,x=1.9,pch=19,col="black",cex=1.5)
points(y=.5958702,x=2.9,pch=19,col="black",cex=1.5)

points(y=.5942029,x=1,pch=15,col="black",cex=1.5)
points(y=.4387755,x=2,pch=15,col="black",cex=1.5)
points(y=.6170213,x=3,pch=15,col="black",cex=1.5)

points(y=.4673913,x=1.1,pch=17,col="black",cex=1.5)
points(y=.3619792,x=2.1,pch=17,col="black",cex=1.5)
points(y=.6136364,x=3.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.4413588,y1=.5390334,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.3468254,y1=.4632105,lwd=2, col="black")
segments(x0=2.9,x1=2.9,y0=.5527818,y1=.6389586,lwd=2, col="black")

segments(x0=1,x1=1,y0=.5458939,y1=.6425119,lwd=2, col="black")
segments(x0=2,x1=2,y0=.3768432,y1=.5007078,lwd=2, col="black")
segments(x0=3,x1=3,y0=.5762737,y1=.6577689,lwd=2, col="black")

segments(x0=1.1,x1=1.1,y0=.4075203,y1=.5272624,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.3096389,y1=.4143194,lwd=2, col="black")
segments(x0=3.1,x1=3.1,y0=.5714921,y1=.6557806,lwd=2, col="black")
legend("bottomright", legend=c("No Rationale","Bureaucratic", "Politicized"), pch=c(19,15,17),
       col=c("black", "black","black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,xlim=c(0,1),ylim=c(0,.8),ylab="Job Performance", col="#0072ce",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy",
     cex.axis=1.5,cex.lab=2)
#axis(1, at=c(0,.2,.4,.6,.8,1),labels=c("0.00","0.20","0.40","0.60","0.80","1.00"),cex.axis=1.5)
lines(x=seq(0,1,by=.05),y=c(0.5606361,
                            0.5556142,
                            0.5505924,
                            0.5455706,
                            0.5405488,
                            0.535527,
                            0.5305052,
                            0.5254833,
                            0.5204615,
                            0.5154397,
                            0.5104179,
                            0.5053961,
                            0.5003742,
                            0.4953524,
                            0.4903306,
                            0.4853088,
                            0.480287,
                            0.4752652,
                            0.4702433,
                            0.4652215,
                            0.4601997),lwd=2,lty=3,col="gray30")
lines(x=seq(0,1,by=.05),y=c(0.5769108,
                            0.579918,
                            0.5829253,
                            0.5859325,
                            0.5889398,
                            0.591947,
                            0.5949543,
                            0.5979615,
                            0.6009688,
                            0.603976,
                            0.6069833,
                            0.6099905,
                            0.6129978,
                            0.616005,
                            0.6190123,
                            0.6220195,
                            0.6250268,
                            0.628034,
                            0.6310413,
                            0.6340485,
                            0.6370558),lwd=2)
lines(x=seq(0,1,by=.05),y=c(0.5922578,
                            0.5757574,
                            0.5592569,
                            0.5427565,
                            0.5262561,
                            0.5097556,
                            0.4932552,
                            0.4767548,
                            0.4602543,
                            0.4437539,
                            0.4272535,
                            0.410753,
                            0.3942526,
                            0.3777522,
                            0.3612517,
                            0.3447513,
                            0.3282508,
                            0.3117504,
                            0.29525,
                            0.2787495,
                            0.2622491),lwd=2,lty=2,col="gray60")
legend("bottomleft", legend=c("Nonpartisan","Copartisan", "Outpartisan"), lty=c(3,1,2),
       col=c("gray30","black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()

###Vote
pdf("Turk2Vote.pdf",width=17,height=5,paper='special') 
par(mar=c(5.1,6, 4.1, 2.1),mfrow=c(1,4))
#Control |    .239521   .0190867      .2020816    .2769603
#Bureaucratic |       .336   .0211448      .2945235    .3774765
#Politicized |   .2111554   .0182338      .1753889    .2469219
plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Rationale",cex.main=2,xaxt="n", xlab="Rationale Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("No Rationale","Bureaucratic","Politicized"),cex.axis=1.5,las=1)
points(y=.239521,x=1,pch=19,col="black",cex=1.5)
points(y=.336,x=2,pch=15,col="black",cex=1.5)
points(y=.2111554,x=3,pch=17,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.2020816,y1=.2769603,lwd=2, col="black")
segments(x0=2,x1=2,y0=.2945235,y1=.3774765,lwd=2, col="black")
segments(x0=3,x1=3,y0=.1753889,y1=.2469219,lwd=2, col="black")

#Control |   .2560706   .0205294      .2157972    .2963441
#Copartisan |    .443418   .0239017       .396529    .4903071
#Outpartisan |   .1555556   .0171043      .1220013    .1891098

plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Partisanship",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
points(y=.2560706,x=1,pch=19,col="black",cex=1.5)
points(y=.1555556,x=2,pch=19,col="black",cex=1.5)
points(y=.443418,x=3,pch=19,col="black",cex=1.5)
segments(x0=1,x1=1,y0=.2157972,y1=.2963441,lwd=2, col="black")
segments(x0=2,x1=2,y0=.1220013,y1=.1891098,lwd=2, col="black")
segments(x0=3,x1=3,y0=.396529,y1=.4903071,lwd=2, col="black")

plot(NA,ylab="Vote for Incumbent?", col="black",
     main="Partisanship by Rationale",cex.main=2,xaxt="n", xlab="Partisanship Treatment", xlim=c(.75,3.25),
     cex.axis=1.5,cex.lab=2,ylim=c(0,.8),pch=19,cex=1.5)
axis(1, at=c(1,2,3),labels=c("Nonpartisan","Outpartisan","Copartisan"),cex.axis=1.5,las=1)
#CC |   .1870968   .0314262      .1254466    .2487469
#CCo |    .427673   .0393595      .3504598    .5048861
#COth |   .1278195   .0290613      .0708088    .1848303
#BC |   .3758389   .0398124      .2977372    .4539406
#BCo |   .5217391   .0426775      .4380169    .6054614
#BOth |   .2101911   .0326216      .1461958    .2741864
#PC |   .2080537    .033366      .1425981    .2735093
#PCO |   .3823529    .041825      .3003031    .4644027
#POth |       .125   .0262277       .073548     .176452

points(y=.1870968,x=.9,pch=19,col="black",cex=1.5)
points(y=.1278195,x=1.9,pch=19,col="black",cex=1.5)
points(y=.427673,x=2.9,pch=19,col="black",cex=1.5)

points(y=.3758389,x=1,pch=15,col="black",cex=1.5)
points(y=.2101911,x=2,pch=15,col="black",cex=1.5)
points(y=.5217391,x=3,pch=15,col="black",cex=1.5)

points(y=.2080537,x=1.1,pch=17,col="black",cex=1.5)
points(y=.125,x=2.1,pch=17,col="black",cex=1.5)
points(y=.3823529,x=3.1,pch=17,col="black",cex=1.5)
segments(x0=.9,x1=.9,y0=.1254466,y1=.2487469,lwd=2, col="black")
segments(x0=1.9,x1=1.9,y0=.0708088,y1=.1848303,lwd=2, col="black")
segments(x0=2.9,x1=2.9,y0=.3504598,y1=.5048861,lwd=2, col="black")

segments(x0=1,x1=1,y0=.2977372,y1=.4539406,lwd=2, col="black")
segments(x0=2,x1=2,y0=.1461958,y1=.2741864,lwd=2, col="black")
segments(x0=3,x1=3,y0=.4380169,y1=.6054614,lwd=2, col="black")

segments(x0=1.1,x1=1.1,y0=.1425981,y1=.2735093,lwd=2, col="black")
segments(x0=2.1,x1=2.1,y0=.073548,y1=.176452,lwd=2, col="black")
segments(x0=3.1,x1=3.1,y0=.3003031,y1=.4644027,lwd=2, col="black")
legend("bottomright", legend=c("No Rationale","Bureaucratic", "Politicized"), pch=c(19,15,17),
       col=c("black", "black","black"), lty=1, lwd=2, cex=1,pt.cex=c(1.25,1.5))

plot(NA,xlim=c(0,1),ylim=c(0,.8),ylab="Vote for Incumbent?", col="#0072ce",
     main="Predicted Values",cex.main=2,xlab="Federal Court Legitimacy",
     cex.axis=1.5,cex.lab=2)
#axis(1, at=c(0,.2,.4,.6,.8,1),labels=c("0.00","0.20","0.40","0.60","0.80","1.00"),cex.axis=1.5)
lines(x=seq(0,1,by=.05),y=c(0.344806,
                            0.3364145,
                            0.328023,
                            0.3196315,
                            0.3112401,
                            0.3028486,
                            0.2944571,
                            0.2860656,
                            0.2776741,
                            0.2692826,
                            0.2608912,
                            0.2524997,
                            0.2441082,
                            0.2357167,
                            0.2273252,
                            0.2189337,
                            0.2105423,
                            0.2021508,
                            0.1937593,
                            0.1853678,
                            0.1769763),lwd=2,lty=3,col="gray30")
lines(x=seq(0,1,by=.05),y=c(0.3555662,
                            0.363037,
                            0.3705077,
                            0.3779785,
                            0.3854493,
                            0.3929201,
                            0.4003909,
                            0.4078617,
                            0.4153325,
                            0.4228033,
                            0.430274,
                            0.4377448,
                            0.4452156,
                            0.4526864,
                            0.4601572,
                            0.467628,
                            0.4750988,
                            0.4825696,
                            0.4900403,
                            0.4975111,
                            0.5049819),lwd=2)
lines(x=seq(0,1,by=.05),y=c(0.2596294,
                            0.2509757,
                            0.2423221,
                            0.2336684,
                            0.2250148,
                            0.2163611,
                            0.2077074,
                            0.1990538,
                            0.1904001,
                            0.1817465,
                            0.1730928,
                            0.1644391,
                            0.1557855,
                            0.1471318,
                            0.1384781,
                            0.1298245,
                            0.1211708,
                            0.1125172,
                            0.1038635,
                            0.0952098,
                            0.0865562),lwd=2,lty=2,col="gray60")
legend("bottomleft", legend=c("Nonpartisan","Copartisan", "Outpartisan"), lty=c(3,1,2),
       col=c("gray30","black", "gray60"), lwd=2, cex=1,pt.cex=c(1.25,1.5))
dev.off()

####REPLICATION FOR APPENDIX D
library(haven);library(interflex)

#####################
##CCES (Figure D1)
#####################
cces_data <- read_dta("CCES18.dta")
ccesModel_Rationale <- lm(dvscale~copartisan*bur*lessIndepend+ideoPre+black+hispanic+educ+gender+age+ownhome+ideo5_01+bornagain+scjob, cces_data)
summary(ccesModel_Rationale)

ccesModel <- lm(dvscale~copartisan*lessIndepend+ideoPre+black+hispanic+educ+gender+age+ownhome+ideo5_01+bornagain+scjob, cces_data)
summary(ccesModel)
d2 <- data.frame(DV=cces_data$dvscale,Treat=cces_data$copartisan,
                 Legitimacy=cces_data$lessIndepend,
                 Rationale=cces_data$bur,
                 ideoPre=cces_data$ideoPre,
                 black=cces_data$black,
                 hispanic=cces_data$hispanic,
                 educ=cces_data$educ,
                 gender=cces_data$gender,
                 age=cces_data$age,
                 ownhome=cces_data$ownhome,
                 ideo5_01=cces_data$ideo5_01,
                 bornagain=cces_data$bornagain,
                 scjob=cces_data$scjob)
out <- interflex(nbins=4,Y = "DV", D = "Treat", X = "Legitimacy", Z = c("ideoPre","black","hispanic","educ","gender","age","ownhome","ideo5_01","bornagain","scjob"), data = d2, estimator = "binning", theme.bw = TRUE,na.rm=T)
pdf("ccesInterflex.pdf",width=6,height=6,paper='special') 
plot(out,theme.bw = TRUE, show.grid = FALSE,Xdistr = "density",
     bin.labs = FALSE,xlim=c(0,1),ylim=c(-.2,.4),
     ylab="Marginal Effect of Copartisan Treatment",
     xlab="Federal Court Legitimacy",main="Copartisan v. Outpartisan",cex.main=1.4)
dev.off()

#####################
##Aug 18 (Figure D2)
#####################
Aug18Turk_data <- read_dta("July18Mturk.dta")

Aug18TurkModel_Rationale <- lm(dvscale~factor(copartisan)*factor(reason)*fjpre+fjaboutright+know+black+hisp+educ+gender+age+ownhome+ideo+fjconf, Aug18Turk_data)
summary(Aug18TurkModel_Rationale)
Aug18TurkModel <- lm(dvscale~factor(copartisan)*fjpre+fjaboutright+know+black+hisp+educ+gender+age+ownhome+ideo+fjconf, Aug18Turk_data)
summary(Aug18TurkModel)
d2 <- data.frame(DV=Aug18Turk_data$dvscale,
                 Treat=Aug18Turk_data$copartisan,
                 Rationale=Aug18Turk_data$reason,
                 Legitimacy=Aug18Turk_data$fjpre,
                 Condition=factor(Aug18Turk_data$allTreat),
                 fjaboutright=Aug18Turk_data$fjaboutright,
                 know=Aug18Turk_data$know,
                 black=Aug18Turk_data$black,
                 hisp=Aug18Turk_data$hisp,
                 educ=Aug18Turk_data$educ,
                 gender=Aug18Turk_data$gender,
                 age=Aug18Turk_data$age,
                 ownhome=Aug18Turk_data$ownhome,
                 ideo=Aug18Turk_data$ideo,
                 fjconf=Aug18Turk_data$fjconf)
d1 <- subset(d2,Treat!=2)
out <- interflex(nbins=4,base="0",Y = "DV", D = "Treat", X = "Legitimacy", Z = c("fjaboutright","know","black","hisp","educ","gender","age","ownhome","ideo","fjconf"), data = d1, estimator = "binning", theme.bw = TRUE,na.rm=T)
pdf("Aug18Turk_CopNon.pdf",width=6,height=6,paper='special') 
plot(out,theme.bw = TRUE, show.grid = FALSE,Xdistr = "density",
     bin.labs = FALSE,xlim=c(0,1),ylim=c(-.2,.4),
     ylab="Marginal Effect of Copartisan Treatment",
     xlab="Federal Court Legitimacy",main="Copartisan v. Nonpartisan",cex.main=1.4)
dev.off()

d1 <- subset(d2,Treat!=1)
out <- interflex(nbins=4,base="0",Y = "DV", D = "Treat", X = "Legitimacy", Z = c("fjaboutright","know","black","hisp","educ","gender","age","ownhome","ideo","fjconf"), data = d1, estimator = "binning", theme.bw = TRUE,na.rm=T)
pdf("Aug18Turk_OutNon.pdf",width=6,height=6,paper='special') 
plot(out,theme.bw = TRUE, show.grid = FALSE,Xdistr = "density",
     bin.labs = FALSE,xlim=c(0,1),ylim=c(-.2,.4),
     ylab="Marginal Effect of Outpartisan Treatment",
     xlab="Federal Court Legitimacy",main="Outpartisan v. Nonpartisan",cex.main=1.4)
dev.off()

d1 <- subset(d2,Treat!=0)
out <- interflex(nbins=4,base="2",Y = "DV", D = "Treat", X = "Legitimacy", Z = c("fjaboutright","know","black","hisp","educ","gender","age","ownhome","ideo","fjconf"), data = d1, estimator = "binning", theme.bw = TRUE,na.rm=T)
pdf("Aug18Turk_CopOut.pdf",width=6,height=6,paper='special') 
plot(out,theme.bw = TRUE, show.grid = FALSE,Xdistr = "density",
     bin.labs = FALSE,xlim=c(0,1),ylim=c(-.2,.4),
     ylab="Marginal Effect of Copartisan Treatment",
     xlab="Federal Court Legitimacy",main="Copartisan v. Outpartisan",cex.main=1.4)
dev.off()

#####################
##Aug 19 (Figure D3)
#####################
Aug19Turk_data <- read_dta("Aug19MTurk.dta")
Aug19TurkModel_Rationale <- lm(dvscale~factor(copartisanCondition)*factor(reason)*fjlegit+female+ideo+know+age+black+hisp+educ+fjjob+fjlib, Aug19Turk_data)
summary(Aug19TurkModel_Rationale)

Aug19TurkModel <- lm(dvscale~factor(copartisanCondition)*fjlegit+female+ideo+know+age+black+hisp+educ+fjjob+fjlib, Aug19Turk_data)
summary(Aug19TurkModel)
d2 <- data.frame(DV=Aug19Turk_data$dvscale,
                 Treat=Aug19Turk_data$copartisanCondition,
                 Rationale=Aug19Turk_data$reason,
                 Legitimacy=Aug19Turk_data$fjlegit,
                 female=Aug19Turk_data$female,
                 ideo=Aug19Turk_data$ideo,
                 know=Aug19Turk_data$know,
                 age=Aug19Turk_data$age,
                 black=Aug19Turk_data$black,
                 hisp=Aug19Turk_data$hisp,
                 educ=Aug19Turk_data$educ,
                 fjjob=Aug19Turk_data$fjjob,
                 fjlib=Aug19Turk_data$fjlib)
d1 <- subset(d2,Treat!=2)
out <- interflex(nbins=4,base="0",Y = "DV", D = "Treat", X = "Legitimacy", Z = c("female","ideo","know","age","black","hisp","educ","fjjob","fjlib"), data = d1, estimator = "binning", theme.bw = TRUE,na.rm=T)
pdf("Aug19Turk_CopNon.pdf",width=6,height=6,paper='special') 
plot(out,theme.bw = TRUE, show.grid = FALSE,Xdistr = "density",
     bin.labs = FALSE,xlim=c(0,1),ylim=c(-.4,.6),
     ylab="Marginal Effect of Copartisan Treatment",
     xlab="Federal Court Legitimacy",main="Copartisan v. Nonpartisan",cex.main=1.4)
dev.off()

d1 <- subset(d2,Treat!=1)
out <- interflex(nbins=4,base="0",Y = "DV", D = "Treat", X = "Legitimacy", Z = c("female","ideo","know","age","black","hisp","educ","fjjob","fjlib"), data = d1, estimator = "binning", theme.bw = TRUE,na.rm=T)
pdf("Aug19Turk_OutNon.pdf",width=6,height=6,paper='special') 
plot(out,theme.bw = TRUE, show.grid = FALSE,Xdistr = "density",
     bin.labs = FALSE,xlim=c(0,1),ylim=c(-.4,.6),
     ylab="Marginal Effect of Outpartisan Treatment",
     xlab="Federal Court Legitimacy",main="Outpartisan v. Nonpartisan",cex.main=1.4)
dev.off()

d1 <- subset(d2,Treat!=0)
out <- interflex(nbins=4,base="2",Y = "DV", D = "Treat", X = "Legitimacy", Z = c("female","ideo","know","age","black","hisp","educ","fjjob","fjlib"), data = d1, estimator = "binning", theme.bw = TRUE,na.rm=T)
pdf("Aug19Turk_CopOut.pdf",width=6,height=6,paper='special') 
plot(out,theme.bw = TRUE, show.grid = FALSE,Xdistr = "density",
     bin.labs = FALSE,xlim=c(0,1),ylim=c(-.4,.6),
     ylab="Marginal Effect of Copartisan Treatment",
     xlab="Federal Court Legitimacy",main="Copartisan v. Outpartisan",cex.main=1.4)
dev.off()

